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Abstract 

The classic Maxwell formula calculates the length of a planar locally 
minimal binary tree in terms of coordinates of its boundary vertices and 
directions of incoming edges. However, if an extreme tree with a given 
topology and a boundary has degenerate edges, then the classic Maxwell 
formula cannot be applied directly, to calculate the length of the extreme 
tree in this case it is necessary to know which edges are degenerate. In 
this paper we generalize the Maxwell formula to arbitrary extreme trees 
in a Euclidean space of arbitrary dimension. Now to calculate the length 
of such a tree, there is no need to know either what edges are degenerate, 
or the directions of nondegenerate boundary edges. The answer is the 
maximum of some special linear function on the corresponding compact 
convex subset of the Euclidean space coinciding with the intersection of 
some cylinders. 

Introduction 

The present paper is devoted to investigation of extreme trees in Euclidean 
spaces. These trees attract interest since the class of such trees is a natural 
extension of the set of locally minimal trees and shortest trees (Steiner minimal 
trees) . The latter ones can be considered as solutions to Transportation problem 
and so have grate importance for applications. It is known, that the searching of 
a shortest tree spanning a given boundary set is a very time-consuming algorith- 
mic problem (an iVP-hard problem), that gives reason for active investigation 
of heuristic solutions. As such solutions, one can choose as extreme trees, so 
as minimal spanning trees. The latter ones are used often, because there is a 
polynomial algorithm of their construction, and quick realizations of it are well- 
known and widespread. The relative error of this heuristic in the worst possible 
situation is called the Steiner ratio of the ambient metric space, see pQ . 

At present, the Steiner ratio is not known for any Euclidean spaces, starting 
with the two-dimensional plane. Notice that in 60th of the previous century, 
E. N. Gilbert and H. O. Pollak [T] conjectured that the Steiner ratio of the 
Euclidean plane is attained at the vertex set of a regular triangle and is equal to 
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y/3/2. But in spite of many attempts of different authors (see a review in [5]), 
this conjecture is not proved yet. The most known attempt was taken in 90th 
by D. Z. Du and F. K. Hwang [3]. But it turns out that their proof contains 
serious gaps which were pointed out as by the authors, so as by other specialists. 
As a result, currently the conjecture is considered as open. 

For the Euclidean spaces of dimension three and more we know even less. 
It is proved in [4] that in these spaces the Steiner ratio is not achieved at the 
vertex set of a regular simplex. Moreover, a fast growing low estimate for a 
possible number of points in boundary set, the Steiner ratio could be achieved 
at is found. Therefore, there is no any reasonable conjecture concerning the 
Steiner ratio value for these spaces. 

In the present paper a new formula is obtained, which gives an opportunity 
to calculate the length of an extreme network spanning a given boundary set 
without the network construction. It turns out that the length of such network 
can be found as a maximal value of some linear function p on an appropriate 
convex compact subset S of the configuration space R N . The function p depends 
on the coordinates of the boundary points only, and the subset S is completely 
defined in terms of the structure of the parameterizing tree of the extreme 
network. It seems to us, that the formula obtained gives a new view onto 
relations between the lengths of distinct extreme trees, in particular, the lengths 
of shortest trees and minimal spanning trees. 

The authors like to use the opportunity to express their gratitude to aca- 
demician A. T. Fomenko for his permanent attention to their work. 

The work is partly supported by RFBR (project N 10-01-00748), President 
Program "Leading Scientific Schools of RF" (project NSh-3224.2010.1), Euler 
program of DA AD, and also by Federal Programs RNP 2.1.1.3704, FCP 02.740.11.5213 
and FCP 14.740.11.0794. 

1 Preliminaries 

Consider an arbitrary tree G = (V,E) with the vertex set V = {vk} containing 
a fixed subset B = {v\, . . . , v n }, and with the edge set E. We assume that all 
the vertices of the tree G having degree 1 or 2 lie in B. Such a set B is called 
a boundary of the tree G and is denoted by dG. Vertices from dG and edges 
incident to such vertices are said to be boundary, and the remaining vertices 
v n+ i, . . . , v n+s and the remaining edges are said to be interior. A network of 
type G in the Euclidean space R m is an arbitrary mapping T V — > R m . The 
restriction of T onto dG is called the boundary of the network T and is denoted 

by ar. 

Each network T is useful to represent as a linear graph, associating each edge 
VkVi of the tree G with the segment [r(ufc),r(wj)l (which could be degenerate). 
This segment is called by edge of the network T. Thereby, an angle between 
adjacent non-degenerate edges of the network is naturally defined. The length of 
the network is also naturally defined as the sum of the lengths of all its edges. A 
network, all whose edges are non-degenerate, is called non- degenerate. Besides, 
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an edge of a network is called boundary (interior), if such is the corresponding 
edge of the tree G. An edge of the tree G is said to be Y -degenerate (T-non- 
degenerate) , if such is the corresponding edge of the network V. 

Let a = {Gk} be a family of non- intersecting subtrees of G. Define the 
reduced tree G/a as follows. Its vertices are all the Gk together with the vertices 
from G, which are not in Gk- Join the vertices from G/a by an edge, if and 
only if they are joined in G by an edge. The boundary of the tree G/a consists 
of all its vertices containing elements from dG. 

For each network T, define the T -reduction of the tree G taking the Gk to 
be equal to the connected components of the subgraph of G generated by T- 
degenerate edges of the tree G. Notice that the network T generates naturally 
the non-degenerate network of the type G/a, which is denoted by V and called 
reduced. 

Let ip B — ► R m be a mapping, which is one-to-one with its image. Consider 
all possible networks with the boundary (p. Then a network having the least 
possible length among all such networks is called a Steiner minimal tree or a 
shortest network with the given boundary. 

A network T is said to be locally minimal, if dT is one-to-one with its image, 
and the angles between adjacent edges of the reduced network T are at least 
120°. Notice that the vertex degrees of the network T can be equal to 1, 2, or 3, 
and all its vertices of degree 1 and 2 are boundary. A locally minimal network 
whose set of boundary vertices coincides with the set of its vertices of degree 1 
is called binary. 

Each shortest network is locally minimal. Locally minimal networks, in 
turn, are shortest "in small", i.e., any sufficiently small part of such network is a 
shortest network with the corresponding boundary. Full information concerning 
locally minimal networks can be found in [5] and [5|- 

For an arbitrary network T of type G we put Zk — Y(vk) £ R™. Let us 
identify the boundary of the network T with the point z = (zi, . . . ,z n ) £ R mn 
of the configuration space R mn . 

By (•,•) we denote the Euclidean scalar product in R m ™. To start with, 
let n = 2, z\ ^ Z2, and v be the direction of the segment [21,22] oriented to 
zi (correspondent, —v is the direction of this segment oriented to z\). In the 
configuration space we define the vector 9 = (—v,v). That is, the vector 9 is 
composed from the directions of the segment [z\, z^\ oriented to its consecutive 
vertices. Then 



i.e. (z, 9) is equal to the length of the segment \z\, Z2]. 

Now, let r be a non-degenerate locally minimal binary network with a bound- 
ary ip. By p(T) we denote the length of the network T. Let Vk be the direction 
of the boundary edge of the network T entering into the point Zk, k < n. Let us 
define the vector 9 of the boundary edges directions of the network T as follows: 
9 = (v\, . . . , v n ). Summing the above expressions for the length of a segment 
and taking into account the additivity of the scalar product and the fact, that 
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the sum of the direction vectors of the edges entering into a vertex of degree 3 
is equal to zero, we have: 

P(T) = (M>. 

The latter equality is referred as Maxwell formula (see, for example, [1], [2], 
|6J). Maxwell formula can be naturally generalized to locally minimal networks 
of general form. In this case the corresponding reduced network can contain 
boundary vertices of degree 2 or 3, where the edges meet by angles of at least 
120°. In Maxwell formula for such networks, at each boundary vertex we need 
to take the sum of vectors of all the entering edges. 

Notice that for non-degenerate locally minimal binary networks in the plane, 
the vector 9 can be calculated by the direction of a single edge and by the 
planar structure of the network. It is convenient to do this in terms of so- 
called twisting numbers. Let us recall the corresponding definition. Let the 
standard orientation of the plane be fixed, that gives us a possibility to define 
the positive (left) and negative (right) rotations. Let e\ and e 2 be two edges of 
a planar immersed binary tree G", and 7 be the unique path in G joining them. 
Then the twisting number tw(ei,e2) from the edge e\ to the edge e 2 is defined 
as the difference between the numbers of left (positive) and right (negative) 
turns in the interior vertices of the path 7 during the motion from e\ to e 2 in 
G. Notice that this function is skew-symmetric and additive along the paths 
(see or 0). 

Let us identify the plane R 2 with the complex field C: the points Zk are 
considered as complex numbers, and the directions v are considered as unit 
complex numbers e 1 ^ . Then the configuration space can be identified with C™ . 
Instead of the Euclidean scalar product, here we consider the Hermitian scalar 
product which is denoted in the same way. 

Let ek be the unique edge of the network T entering into the point Zk, 
k < n, and 9 = (e 1 ^ 1 , . . . , e 1 ^™) be the directions vector of the edges e^. Then 
the number (z, 9) is real, and it is equal to the length p(T) of the network T. 
Indeed, it suffices to verify this equality for a network having a single edge 
[zi,z 2 ], z\ 7^ z 2 . Let e 1 ^ be the direction of the segment [zi,z 2 ] oriented from 
z\ to z 2 (in this case — e 1 ^ is the direction of this segment oriented to the point 
zi). Then 9 = (-e^,e^), and hence 

(z, 9) = zi(-e^)+z 2 e^ = (z 2 -zi)e~ 1 ^ = {z 2 -zi)(z 2 - Zi)/\z2-z\\ = \z 2 -zi\. 

Let t pq = tw(e p , e q ) be the twisting number from the edge e p to the edge e q . 
Then 

Put tk = (e 1 3"* fcl , . . . , e* 3" tfc ") . Notice, that the fcth component of the vector 
ife is equal to 1. Then 9 — —e l ^ k tk- Thus, starting with the direction of one 
boundary edge only, one can find the directions of all the remaining boundary 
edges (and all the other edges also) of the tree T using just the twisting numbers. 

As an application of the above formula, let us calculate the length of the 
network T and the directions of all its edges without an explicit construction 
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of the tree and using only the information concerning the boundary mapping 
tp and the planar structure of the network T (the twisting numbers for all the 
pairs of edges of the corresponding immersed planar tree G'). Since 

p{T) = {z,6) = -e i ^{z,t k ), 

then p(T) = \(z,tk)\ and ipk is equal to the argument of the number — {z, 6)/{z,tk), 
and since {z, 0) is a positive real, then ipf. is equal to the argument of the number 
— (tk, z). If we know ipf., then we can find out the remaining ip p by means of the 
above formulas. 

To calculate similarly the length and the edges directions of a locally mini- 
mal planar network of general form, we need to partition the reduced network 
r into the binary components, cutting it by the vertices of degree 2 and bound- 
ary vertices of degree 3, and to proceed the calculation for each component 
separately 

Let us return to the networks in W m . By [G,tp] we denote the set of all 
the networks in R m having the type G and the boundary ip. Each network T 
from [G, ip] is uniquely defined by the images Zk = T(yk), k > n, of all the 
interior vertices of G. Consequently writing down the vectors z k , k > n, as the 
components of the vector (z„ + i, . . . , z n+s ), we identify the set [G, p] with the 
space R ms . Notice, that the length of the network is a real- valued function pg, v 
on [G, ip] . It is easy to see, that this function is convex and tends to infinity 
as the arguments unlimitedly increase. Therefore, the set of minima of this 
function is non-empty and convex. Each network corresponding to a minimum 
of this function is said to be an extreme network of the type G with the boundary 

<P- 

Recall that if among the extreme networks of a type G with a boundary tp 
a locally minimal network exists, then the extreme network in [G, tp] is unique, 
see [2]. Thus, [G, tp] contains at most one locally minimal network. Notice, that 
an extreme network need not always be locally minimal. For example, it can 
contain vertices of degree more than 3. 

By definition, the length of an extreme network can be calculated as the least 
value of the function pc,<p on ffi. ms . But this function is a sum of square roots of 
square polynomials on the coordinates of the interior vertices, that complicates 
the investigation of the extreme networks in its terms. In the present paper we 
show, how the length of an extreme network can be calculated by maximization 
of a linear function. But this maximization need to be proceeded on a more 
complicated subset of the configuration space, namely, on an intersection of some 
cylinders and a linear subspace. We generalize the Maxwell formula and make it 
uniform for all types of networks obtained from a given type G by degeneration 
of some edges of the graph G. Besides, in the case of planar networks, we escape 
the necessity to go over all possible non-equivalent immersions G' of a binary 
tree G. 
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2 Generalized Maxwell formula 

Let G = (V, E) be an arbitrary tree with a boundary B = {v\, . . . ,v n } C V 
consisting of n elements. Using the structure of the tree G, let us form a system 
Sq consisting of equations and inequalities on variables 9 3 k which are considered 
as the standard coordinates (6{, . . . , 0f , . . . , 0* , . . . , 6™) in the space R mn . We 

put e k = {el,..., e™) ande = {e 1 ,...,e n ). 

Let e be some edge of the tree G. By G r = (V r ,E r ), r = 1, 2, we denote 
the connected components of the graph G\e. Thus, G \ e = G\ U G2, and put 
B r = B n V r . By Voie) we denote the resulting partition {Bi,^} of the set 
B. Let us choose one of B r £ Vc(e), and let B r — {v^, . . . ,Wfe p }- By o~ e we 



_ 2 

denote the inequality J2q=i @k„ < L and by cr we denote the vector equation 

J2k=i ®k = 0. The system Sq is formed from a and o e over al edges e of the 
tree G. 

Let \Sg\ C R" m be the set of all solutions to the system Sq- Notice that 
\Sg\ does not depend on the choice of the components Bk, since equality is 
valid. Besides, each inequality a e defines a convex subset in R mn , bounded by 
an elliptic cylinder, i.e. it is the product of an m-dimensional elliptic disk and 
l>m(n-i). The origin together with the points of the subspace II defined by 
equation <r, which are close to 0, are solutions to the system Sg- Therefore, 
\Sq\ is a convex body in the subspace II. 

Let fB — > R m be an arbitrary mapping, and z k — (p(vk), and z = (z\, ... , z n ). 
Put p v {0) = (6,z). 

Theorem 1 Under the above notations, the length of each extreme network 
from [G, ip] is equal to the largest value of the linear function p v on the convex 
set\S G \. 

Proof. Let 8 £ \Sg\ be an arbitrary point of the set \Sq\- To each pair (v,e), 
where e £ E is the edge of the tree G incident to the vertex v £ V, we assign 
the vector 9(v,e) £ R m using the vector 9 as follows. Let Pg(^) = {-817^2}) 
and v £ Bi and B\ = {vk t , • • • , Vk p }■ Then we put 

p 

9=1 

Lemma 1 For any edge e = vw of the tree G we have 9(v, e) — —9(w, e). 
Proof. The statement of Lemma follows immediately from equality a. 
Lemma 2 For any boundary vertex Vk, 1 < k < n, we have 9k — 0(vk,e). 

e:v k £e 

Proof. Let e q = VkW q , q = l,...,p, be all the edges from G, which are incident 
to Vk- We cut the tree G by the vertex Vk, and let G q = (V qi E q ) be the 
component containing w q . Put B q = V q n B, then 

B = B 1 U • • • U B p U {v k }, 
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hence 

p p 

= J2 + °k = - 9 ( V k, e ?) + Ok, 

q=l q=l 

where the last equality follows from Lemma [TJ which was to be proved. 

Lemma 3 Let ej, j = 1, . . . , r, be all the edges of the tree G, which are incident 
to its interior vertex v. Then the equality 

r 

holds. 

Proof. Indeed, let ej = vvjj, then, due to Lemma [TJ we have: 9(v,ej) = 
-9(wj,ej). Let V G (ej) = {B{,B J 2 }, and v G B{ for each j. Then B = UjB J 2l 
and therefore, 

n 

j j k=l 

which was to be proved. 

Lemma 4 For each vertex v of the tree G and each its edge e incident to v, we 
have \\0(v,e)\\ < 1. 

Proof. Taking into account equality a (9) , the statement of Lemma is equivalent 
to the validity of inequality cr e {9). 

Let L £ [G, if] e an extreme network, and {i>„+i, . . . , i> n + s } be the set of all 
interior vertices of the tree G. For n + l<i<n + swe also define Zk as follows: 

Zk = r(wfc). 

Lemma 5 Under the above notations, we have 

(z,9) = ]T (z k ,9(v k ,e)). 

v k Ge 

Proof. Let us partition the sum in the right hand part of the equality into two 
sums: the first one is over all interior vertices v k , and the second one is over all 
the boundary vertices. The first sum vanishes due to Lemma [3J In the second 
sum we group the terms corresponding to the same vertex and apply Lemma [2j 
Lemma is proved. 

Thus, due to Lemma [SJ 

(z,9)= Y ( z k,0(vk,e)) = Y ( z k~ zi,9(v k ,e)), 

{v k ,e) e=v k v l 
v k &e 
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therefore, due to Lemma 01 we have 

{z,0) = E (zk -zi,9(v k ,e)) < E \\ z k- zi\\= p{Y). 

e=v k v t e=v k v t 

Since, as we remember, 9 is an arbitrary point from \Sg\, we conclude that the 
maximal value of the function p v (9) on \Sq\ does not exceed p(T). 

Now, let us show that this maximal value is reached. For each T-nondegenerate 
edge e = vw, by £,(v,e) we denote the unit vector from E m having the same 
direction as the edge e of the network T entering the point T(v) has. Notice that 
£(v, e) = — e). Further, due to Extreme Networks Local Structure Theorem 
(see [2], Theorem 4.1, or [5], Theorem 4.1), for any T-degenerate edge e = vw 
the pair (v,e) can be assigned with some vector £(u,e) € M m , \\£{v, e)|| < 1, in 
such a way, that the equality £(i>, e) = — £(u>, e) holds and also the vector- valued 
function £ (now £ is defined on all the pairs (v, e) where vertex t> is incident to 
edge e) meets relation 

at any interior vertex u of the tree G. We put = X^ e -u fc ee e), l<fc<n 
and show that the vector £ = . . . ,£ n ) g M™ 1 " is a solution to the system Sg 
and^(e) = (z,£)=p(r). _ 

Indeed, consider an arbitrary edge e = vw of the tree G, and let G \ e = 
Gi U G 2 , G r = (F r , # r ), and v € Vi and B x = B n F x = ,...,%,}. 

Lemma 6 Under the above notations, we have £,(v,e) = X)g=i£fc<j- 

Proof. At first, assume that uG-B and u = Ufc s . Since the equality £(v',e') = 
— £(w/, e') holds for any edge e' = i/u>', we conclude that 

°= E (^',e')+eK,e')) = 

= E E £K> e ')+ E ^ e ')+ 

gG{l,...,p}\{s} e' eE 1 :v kq Ee' e'eEi:vee' 

+ E E ^V) = 

d'SVi\-Bi e'GEiiu'ee' 

E &« + (&.-£(»,«o) + E E ^'> e ')- 

ge{i,...,p}\{s} ti'eVi\Bi e'eBnu'ee' 

Since the latter sum vanishes in accordance to equality ([1]) , we obtain the state- 
ment of Lemma for the boundary vertex v. 
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Now, let v be an interior vertex of the tree G. Then 
°= E (£(«', e')+£(«/,eO) = 

= E E £K< e ')+ E ^ e ')+ 

9=1 e':i)fc ge' e'G-Ei:t)Ge' 



+ E E 

« evi\( Biu{«} I 



The first sum in the latter formula coincides with £k q due to definitions, 

the second one is equal to — e) in accordance with equality ([T]), and the third 
one vanishes due to equality (pj also, thus, we get the statement of Lemma for 
the interior vertex v. Lemma is proved. 

Lemma 7 The equality X)L-=i Cfc = holds. 

Proof. Again, taking into account the fact that the equality £ (v, e) = — £,(w, e) 
is valid for any edge e = vw, we get 

n 

0= (av,e)+z(™,e)) = j2 E £K> e )+ E E ^> e )- 

e=vw£E fc=l e:« fc £e t)6(V\i3) e£-E:u£e 

The first sum in the latter equality coincides with J\=i & due to definitions, 
and the second one vanishes due to equality (JlJ. Lemma is proved. 

Lemma|n]and the condition e)|| < 1 imply that £ meets each inequality 
cr e . Lemma [7] implies that £ also meets condition a. Thus, £ € |Sg|. Further, 

n n 

p<p(o = ( z >o = ^2( z k,tk) = E E (^)^( x fe' e )) = 

fe=l k=l e:t) fc 6e 

n+s 

= E E ( z k,£( x k,e)) = E ( Zk ~ z hi( v k,e)) = 

k=l e:v k £e e=v k vt£E 



= E ( Z k-Zh€(vk,e)) = 

e=v k vieE e=VkVi£E 



Zk - Zl 

Zk - zi, ■ 



\\ z k - zi\\ 



= J2 \\zk-z l \\=p{T). 

e—v^vi £E 



Theorem Q] is proved. 



Remark. Notice that, due to Lemmas [5] and 21 each component of vector 9 € 
\Sq\ is bounded, therefore \Sc\ is a compact subset of M. mn . Thus, taking the 
above into account, we conclude that \Sg\ is a convex compact. 
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As above, let G — (V, E) be a tree with a boundary B = {v\, . . . , v n } C 
and Ed C E be some family of edges of the tree G. By a — {Gk} we denote 
the family of nonintersecting subtree of G, such that the union of all their edges 
coincides with Ed- Assume that each Gk intersects B by at most one vertex. 
Such families Ed are said to be admissible. Consider the reduced tree G/a and 
choose its boundary to be the set of all the vertices which intersect the set B. 
Due to the restrictions imposed on the subtree Gk from a, the boundary set 
obtained consists of the same number of points as the set B does, that gives us 
an opportunity to identify it with B. Thus, for an arbitrary mapping tpB — > R m , 
the next two spaces are defined: [G, <p] and [G/a, <p]. If T € [G,<p] degenerates 
all the edges from Ed, then T generates naturally the network from [G/a, p), 
which we denote by T/a. 

As above, consider the system Sg consisting of equation a and inequali- 
ties <7 e , and through out of it all the inequalities corresponding to the edges 
e € Ed- By Sg\ Ed we denote the resulting system. 

Theorem 2 Under the above notations, the length of extreme network T/a G 
[G/a, ip] is equal to the largest value of the linear function p v on the set \ Sc\Ed\ 
of all the solutions to the system Sg \ Ed ■ 

Proof. It is sufficient to notice that Sc\Ed — Sc/ a , since for any edge e from 
G/a we have: Vciaifi) = T-'g(s)- Theorem[2jis proved. 

Thus, throwing out of inequalities of the form a e from the system Sg is 
equivalent to factorization of the tree G by the corresponding edges e. The 
authors hope that the formula obtained for the calculation of extreme networks 
lengths can be useful in investigation of Steiner ratio of Euclidean spaces. 

Remark. The results of this paper could be generalized in the following di- 
rections. First, Extreme Networks Local Structure Theorem from (2] and [5] 
is proved for the case of extreme weighted networks. The classical Maxwell 
formula also can be easily transformed to this case (namely, the length of the 
direction vector of a weighted edge should be equal to the weight of the edge). 
Therefore, the above results can be easily generalized to the case of weighted 
trees. 

Second, the classical Maxwell formula remains valid for extreme networks 
with cycles. Therefore it seems not very difficult to generalize Theorems [1] 
and [2J to the case of arbitrary extreme networks in R m (not necessary the trees 
as in the present paper), and also to the case of arbitrary weighted extreme 
networks in M. m . 

At last, in [2j we obtained theorems on the local structure of extreme net- 
works in normed spaces. It would be interesting to work out an analogue to the 
Maxwell formula for this case in terms of so-called p- impulse, see [2], and after 
that to generalize the results of the present paper to the case of normed spaces. 
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